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Abstract 

Strong wind-wind collisions in massive binaries generate a very hot plasma that frequently produces a moderately strong 
iron line. The morphology of this line depends upon the properties of the wind interaction zone and its orientation with 
respect to the line of sight. As the binary components revolve around their common centre of mass, the line profiles 
are thus expected to vary. With the advent of the next generation of X-ray observatories {Astro-H, Athena) that will 
offer high-resolution spectroscopy above 6keV, it will become possible to exploit these changes as the most sensitive 
probe of the inner parts of the colliding wind interaction. Using a simple prescription of the wind-wind interaction in 
an early-type binary, we have generated synthetic line profiles for a number of configurations and orbital phases. These 
profiles can help constrain the properties of the stellar winds in such binary systems. 
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1. Introduction 

Chandra and XMM-Newton have opened up the era of 
high-resolution X-ray spectroscopy for moderately bright 
X-ray sources. In the case of massive stars, HETG and 
RGS spectra unveiled for the first time some details of the 
morphology of the spectral lines (see e.g. Giidel & Naze 
, 2009). For presumably single massive stars, where the 
X-ray emission arises from hot plasma embedded in the 
stellar wind, a number of observational and theoretical 
studies were undertaken with the goal to connect the line 
morphologies to the properties of the stellar wind (e.g. 
Feldmeier et ah, 2003; Owocki & Cohen, 2006; Oskinova 
et ah, 2007; Cohen et ah, 2010; Herve at ah, 2013; Rauw 
et ah, 2015). 

In massive binary systems, the X-ray emission not only 
arises from the individual winds of both stars, but also 
from shock-heated plasma inside the wind interaction zone 
between the stars (e.g. Stevens et ah, 1992). At the shock 
front between the winds, the kinetic energy normal to the 
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shock is converted into heat, thereby generating a very hot 
plasma. Since the plasma temperature exceeds the values 
typically reached in the shocks that prevail in the winds 
of most single massive stars, the detection of such a hot 
plasma provides strong hints of a wind-wind interaction. 
The most obvious signature of such a high-temperature 
plasma is the moderately strong iron line near 6.7 keV, 
commonly called the Fe K line, that can be observed in 
some systems. Whilst numerous studies have been devoted 
to colliding winds over the last two decades (for a general 
review on colliding winds see Rauw & Naze, 2015), most 
of them focused on broadband medium-resolution spec¬ 
troscopy. Based on density and velocity maps obtained 
from hydrodynamical simulations, Henley et al. (2003) 
presented calculations of theoretical X-ray line profiles of 
colliding wind binaries and of their variations along the 
orbit for a range of orbital and wind parameters. Their 
study focused on the Lya transitions of O viii, Ne X, Mg 
XII, Si XIV and S xvi. These lines are located at energies 
below 3 keV, in the sensitivity band of the current high- 
resolution spectrographs. However, most of these lines are 
also emitted by the plasma in the winds of each individual 
binary component outside the wind-wind interaction zone. 
Their profiles are thus often contaminated, or even domi¬ 
nated, by contributions from the intrinsic emissions of the 
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stars. Moreover, as shown by Henley et al. (2003), photo¬ 
electric absorption by the cool unshocked winds affects the 
morphology of the lines at longer wavelengths, especially 
in the case of the O, Ne, and Mg La lines. Both effects 
render the interpretation of the observed line profiles in 
terms of colliding winds rather difficult. 

The Fe K feature offers a promising way out of this 
dilemma. Indeed, except for some stars featuring a strong 
magnetic field, this line is not seen in the X-ray spectra 
of single massive stars. Furthermore, at energies near 6.7 
- 7.0 keV, it is essentially unaffected by absorption in the 
unshocked winds. Moreover, it is rather isolated, thereby 
limiting the problems of blends with other species. Un¬ 
fortunately though, the iron line falls outside the sensi¬ 
tivity range of the current generation of high-resolution 
spectrographs. However, this situation is about to change 
in the near future thanks to the JAXA mission Astro-H 
(Takahashi et ah, 2014) and, in the more distant future, 
thanks to the ESA observatory Athena (Nandra et ah, 
2013). Each of these missions will carry a high-resolution 
bolometric spectrograph optimized for a high-sensitivity 
and high-resolution coverage of the spectral region around 
6.0 - 7.0 keV. Astro-H will be equiped with the Soft X-ray 
Spectrometer (SXS, Mitsuda et ah, 2014), whilst Athena 
will host the X-ray Integral Field Unit (XIEU, Barret et 
ah, 2013; Ravera et ah, 2014). Whilst these instruments 
are primarily designed to address other scientific questions, 
they can also be efficiently used for the purpose of the 
study of colliding wind binaries (Sciortino et ah, 2013). 

In this context, we address here the question of the 
morphology of the Fe line at 6.7 keV and its variations 
with orbital phase. 

2. The iron line in the spectra of colliding wind 

binaries 

Before we turn to the discussion of our models, we need 
to briefly consider the nature of the Fe K line that is ob¬ 
served in the spectra of colliding wind binaries at ener¬ 
gies around 6.7keV. As shown by spectral fits, this line 
is compatible with emission from highly ionized iron in 
a thermal plasma of high temperature T. The emissiv- 
ity of the triplet of helium-like Fe xxv between 6.64 and 
6.70 keV peaks near logT ^ 7.8 (i.e. 63 MK or kT = 
5.4keV). Conversely, the emissivity of the Lya doublet 
of hydrogen-like Fexxvi peaks at logT ~ 8.2 (i.e. 158 MK 
or kT = 13.7keV). Considering that typical wind veloci¬ 
ties of massive stars range between 1000 and 3000 kms“^, 
the Rankine-Hugoniot condition for strong shocks implies 
typical post-shock temperatures between 1.2 and 10.5 keV. 
This leads us to conclude that the iron line of the major¬ 
ity of the colliding wind binaries should be associated with 
the Fexxv helium-like triplet. 

Helium-like triplets consist of a resonance line, an inter¬ 
combination doublet and a forbidden line. In the case of Fe 
xxv, the energies of these lines, taken from the AtomDB 
database (Foster et al., 2012), are 6.7004keV (resonance 
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Figure 1: Modification of the 'IZ ratio between the forbidden and 
intercombination lines in the Fe xxv triplet formed in the wind of an 
O-type star with Tgff = 55 000 K and log g = 4.0. 


line), 6.6823 and 6.6676keV (intercombination doublet), 
and 6.6366keV (forbidden line). 

As shown e.g. by Blumenthal et al. (1972) and Porquet 
et al. (2001), the forbidden line of the helium-like triplets 
of light elements can be strongly suppressed as a result 
of either collisional excitation or radiative pumping of the 
electrons from its upper level (2^5'i) to the upper levels 
of the intercombination doublet (2^Pi_2)- The latter ef¬ 
fect has been used extensively to study the location of the 
X-ray emitting plasma in the winds of massive stars (e.g. 
Leutenegger et ah, 2006). Does a similar effect also exist 
for the Fexxv triplet? To address this question, we fol¬ 
low the formalism of Blumenthal et al. (1972); Porquet et 
al. (2001) and Leutenegger et al. (2006). These authors 
have shown that the ratio TZ between the forbidden and 
intercombination lines can be written 
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Here TZo is the ratio between the forbidden and intercom¬ 
bination line in the absence of photospheric light and col¬ 
lisional excitation. Ne is the local electron density and 
Nc is the critical electron density above which collisional 
excitation from 2^51 to 2^Pi 2 becomes important. The 
photoexcitation rate due to the photospheric radiation 
field is multiplied by the dilution factor 
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at a distance r from the centre of the star of radius i?*. 
This value is compared to a critical rate ipc- In the case of 
the Fexxv triplet, Blumenthal et al. (1972) quote values 
ofNe = 4.7 X 10i®cm-3 and ife = 1-23 x 10®s-h 
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For a solar composition, Nc corresponds to pc = 2.3 x 
10“®gcm“^. This value is several orders above the typ¬ 
ical density in adiabatic post-shock regions (e.g. Pittard, 
2009). The suppression of the forbidden line by collisional 
excitation can thus be neglected. 

Concerning the impact of photoexcitation, the relevant 
transitions that will pump the electrons are located in 
the EUV at wavelengths of 271.4 and 396.7 A (Kramida 
et ah, 2014). We use the fluxes at the stellar surface 
predicted by the non-LTE plane-parallel TLUSTY stellar 
atmosphere model and available via the Ostar2002 grid 
(Lanz & Hubeny, 2003) to estimate ip*- If we consider a 
very hot main-sequence star with T^s = 55 000 K, we find 
that the impact of photospheric radiation on the Fexxv 
triplet is rather marginal (of order 8% at the stellar sur¬ 
face, see Fig. 1) and drops rapidly when we move out into 
the stellar wind, at radii > 5 7?*, where the wind-wind 
collision occurs. The main reason for this weak effect is 
of course the high value of ipc- even such a hot star does 
not provide enough radiation at the relevant wavelengths 
to produce a strong pumping. We can thus conclude that 
we can safely neglect the impact of photospheric radiation 
on the relative strengths of the components of the Fe xxv 
triplet. 

3. A simplified model for the wind-wind interac¬ 
tion 

Depending on the efficiency of radiative cooling in the 
post-shock plasma, colliding wind binaries are frequently 
divided into two broad categories (Stevens et ah, 1992). 

The plasma is said to be in the adiabatic regime when 
the post-shock densities are sufficiently low for the cooling 
time (icooi) to be much longer than the escape time (tesc) 
from the shock region. The efficiency of radiative cooling is 
thus expressed via the parameter y = tcooi/Asc- This sit¬ 
uation mainly occurs in wide, long-period binaries and/or 
systems consisting of stars with relatively low mass-loss 
rates. The corresponding wind interaction zone is quite 
smooth and can usually be well modelled with hydrody- 
namical simulations. The X-ray emission should scale as 
n-3-^ d~^ where M, v and d are the mass-loss rate, the 
pre-shock wind velocity and the separation between the 
stars, respectively. 

Conversely, if the post-shock plasma density is high, ra¬ 
diative cooling happens very quickly and the wind interac¬ 
tion zone is then said to be radiative. This concerns close 
binary systems hosting stars with large mass-loss rates, 
but can also occur in wide, eccentric systems around peri- 
astron passage. From first principles, the X-ray emission 
of these systems should scale with the kinetic power of 
the incoming stellar winds. However, the observed level of 
X-ray emission is usually much lower than theoretically ex¬ 
pected. This could be due to radiative inhibition (Stevens 
& Pollock, 1994), radiative braking (Cayley et ah, 1997) 
or the strong hydrodynamical instabilities that strongly 
distort the wind interaction zone (Kee et ah, 2014). 


In our present work, we consider mainly systems that 
can be described as undergoing an adiabatic wind inter¬ 
action. We further assume that the winds collide at their 
terminal velocity Uoo, thereby neglecting wind accelera¬ 
tion, radiative inhibition and braking. Since the emission 
of interest here arises in the inner parts of the interaction 
zone, we neglect the impact of the orbital motion and the 
shocks are considered to be axisymmetric about the bi¬ 
nary axis. We then use the analytical solution of Canto 
et al. (1996) to describe the properties of the wind in¬ 
teraction zone. In principle, the formalism of Canto et 
al. (1996) holds for situations where the shocked plasma 
forms a thin shell around the contact discontinuity, which 
is the case if both winds are highly radiative. Although 
we consider here mainly situations where the shocks are 
adiabatic, the position of the contact discontinuity is un¬ 
likely to drastically differ from the predictions of Canto et 
al. (1996), as was shown by Parkin & Pittard (2008) and 
Pittard (2009). 

For a given value of the wind momentum ratio 

Ml Voo,2 
Ml Voc,l 

solving equation (28) of Canto et al. (1996) by means of a 
Newton-Raphson scheme allows us to estimate the asymp¬ 
totic opening angle of the shock 9oo as seen from the star 
with the weaker wind. The interval [O,0tx3] is then dis¬ 
cretized into 200 steps 9k- For each angle, equations (23) 
and (24) of Canto et al. (1996) are solved to compute the 
shape of the contact discontinuity, i.e. to establish the rela¬ 
tion between r and 9 (see Fig. 2). Along with the adopted 
mass-loss rates and wind velocities, equations (29) and 
(30) of Canto et al. (1996) further allow us to estimate 
the mass surface density cr of the wind interaction zone 
and the tangential velocity along the contact dicontinuity. 
The resulting shape and velocity vector are then rotated 
about the binary axis to generate an axisymmetric shock 
cone, discretized into 200 x 360 2-D cells, along with the 
associated velocity field. The thickness of these cells is 
computed as the ratio of the surface density tr divided 
by the post-shock density ps- Since we are dealing with 
an adiabatic wind interaction zone, we use the Rankine- 
Hugoniot condition for strong shocks ps = ^Pw where pw 
is the pre-shock wind density. The volume associated with 
a cell is then calculated following 

„, sin 9 d9d(p a 

dV = — -^ — (4) 

sm (a — 9) Ps 

(see Fig. 2). Here d(p = 1° is the incremental step of the az¬ 
imuthal angle used to build the axisymmetric shock cone. 
a is the angle between the tangential velocity and the di¬ 
rection of the binary axis. 

The associated contribution to the line emission is given 

by 

de = QnenHdV (5) 


3 




z A 




Figure 2: Geometry of our model of the wind interaction zone. The 
star with the weaker wind is on the left. The contact discontinuity 
between the winds is shown by the thick black line. The blue vector 
indicates the tangential velocity at the point of coordinates (r,d). 
/3 = Tt 12 — a + d is the angle between the tangential velocity and the 
normal to the position vector. 

where and nn are the electron and hydrogen number 
densities in the post-shock region and Q is the line emis- 
sivity at the temperature of the post-shock plasma inside 
the cell assuming ionization equilibrium. The post-shock 
temperature is computed assuming that the velocity nor¬ 
mal to the shock is entirely thermalized (see Stevens et ah, 
1992) and that the electrons and ions have equalized their 
temperatures. For wind-wind interactions in wide binary 
systems, these two points, ionization equilibrium and equal 
electron and ion temperatures (T^ = Tjon), are clearly ap¬ 
proximations. Indeed, as pointed out by several authors 
(e.g. Usov, 1992; Zhekov & Skinner, 2000; Zhekov, 2007; 
Pollock et ah, 2005), the shocked plasma is probably out of 
ionization equilibrium and the shocks are likely collision¬ 
less, implying < Tion- However, these approximations 
are frequently used (including in the work of Henley et 
ah, 2003). Moreover, they mainly affect the strengths of 
the simulated lines and should be much less important for 
their profiles. 

For a given pair of orbital inclination i and orbital 
phase $, we determine the direction of the line-of-sight 
and project the tangential velocity of each cell onto this 
direction to obtain the line-of-sight velocity of the gas in 
the cell. If aberration of the shock due to orbital motion 
can be neglected, then the orientation of the line-of-sight 
can be characterized via a single angle 0 with 

cos(0) = sin(f) cos(M) (6) 

Here M = v— where v is the true anomaly and w is 
the longitude of periastron of the primary star, which we 
consider to be the star with the stronger wind (see Fig. 3). 

Although the effect of photoelectric absorption by the 
winds is very low at the energies of interest here, we have 
nevertheless accounted for the optical depth along the line- 
of-sight. Using the ionized wind opacity model of Naze et 
al. (2004), we estimate the opacity of the wind at 6.7keV 





Figure 3: Illustration of the angles used to define ©. The star with 
the weaker wind is the one at the origin of the axes. The x and y 
axes are inside the orbital plane and rotate as the stars move around 
each other on their orbit. 

to be 0.7 cm^ g“^ for solar abundances. This is about a 
factor 100 and 40 lower than the equivalent opacity at the 
energies of the O viii La and Ne x La lines, respectively. 
To compute the optical depth for a given sightline, we first 
check whether it intersects any of the stars and whether 
it crosses the wind of one or two stars. If the cell is oc¬ 
culted by one of the stars, its contribution to the observed 
profile is set to zero. Otherwise, the optical depth along 
the sightline is evaluated accounting for the appropriate 
opacities of each wind that is crossed. 

Thermal broadening is accounted for by distributing 
the flux of the cell over a Gaussian profile centered on 
the line-of-sight velocity. Given the high mass of the Fe 
ions, its effect remains quite moderate though (typically 
ctherm = 86kms“^ for a 50MK plasma). 

Finally the synthetic line profile is obtained as an his¬ 
togram of the contributions of all the cells according to 
their line-of-sight velocity. We assume a constant energy 
resolution of the instrument of 3eV (FWHM) at the rel¬ 
evant energies. The synthetic profiles are computed with 
line-of-sight velocity steps of 50 km s“^, and are then binned 
onto line-of-sight velocity bins of 150kms“^ corresponding 
to a bit more than one resolution element. 
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Table 1: Parameters of the test models for circular binary systems 
with an orbital separation of 100 Rq. All systems feature the same 
primary star, but differ by the properties of the secondary. 


Parameter 

Primary 

I 

Secondary 

H 

HI 

M (lO-^Mgyr-i) 

1.0 

0.98 

0.3 

0.1 

Voo (kms“^) 

2000 

2000 

2000 

2000 

R* (R-o) 

20 

20 

10 

10 

M, (M 0 ) 

30 

30 

30 

30 

V 

- 

0.98 

0.3 

0.1 

X 

> 5.6 

5.7 

13.1 

26.8 

OooH 


90.4 

111.5 

128.8 


4. Results 

In this section, we present the output of our model for 
several sets of parameters, including systems with circular 
or elliptical orbits. As pointed out above, the Fe xxv 
line consists of four components which are closely spaced 
in wavelength and will thus be blended over parts of the 
parameter space. To illustrate some effects, we thus start 
our discussion with the simulation of a simpler case, where 
we consider a single line. 

4.I. Systems with circular orbits 

Let us consider early-type binaries with circular orbits. 
As a first step, we compare the line profiles obtained for 
three different values of the wind momentum ratio 77 and 
seen under five different viewing angles. To ease compar¬ 
ison with the work of Henley et al. (2003), we adopt the 
same values of 0 as done by these authors. 

Because of numerical singularities our code cannot deal 
with the case where rj is strictly equal to one. Case I thus 
rather considers rj = 0.98. Figured illustrates the profiles 
of individual lines as a function of viewing angle for the 
three cases described in Table 1. Figure 5 illustrates the 
expected morphologies of the full Fe xxv complex. 

We start by considering Case I, which corresponds to 
a situation where the wind interaction zone is almost pla¬ 
nar and located midway between the stars. Therefore, for 
conjunction phases (0 = 0 ° for the conjunction with the 
primary in front and 0 = 180° for the opposite situation), 
the velocity vectors of the material flowing out of the wind 
interaction zone are essentially perpendicular to the line of 
sight. This results in very narrow profiles, centered on zero 
velocity. For these configurations, one hence expects to 
resolve the individual components of the Fe xxv complex 
(see Fig. 5). We note that our results slightly underesti¬ 
mate the line width at these viewing angles. In a real 
system, the velocity vectors in the wind interaction zone 
are not simply tangent to the contact discontinuity, but 
their direction spans a range of values around the tangent. 
The hydrodynamic simulations of Henley et al. (2003) in¬ 
deed yield somewhat wider lines for 77 = 1 and 0 = 0 or 
180°. For other viewing angles. Case I shows very large 
variations of the width of the line. The width is maximum 


for 0 = 90°, i.e. when the line of sight crosses the plane of 
the wind interaction region. For this orientation, the full- 
width of the line could reach 2 x Voo if ah cells of the wind 
interaction zone were emitting at the same level. The fact 
that the line width does not reach this value reflects the 
variations of the plasma temperature and hence emissivity 
along the shock: the further we move away from the binary 
axis, the more oblique the shock with respect to the inflow¬ 
ing wind and thus the lower the post-shock temperature 
and hence the lower the contribution to the line emission. 
We further see that the profile remains roughly symmetric 
and centered on zero velocity. Unlike the lines simulated 
by Henley et al. (2003), the profiles in Fig. 4 do not dis¬ 
play strong asymmetries between the blue and red side. 
This is because the impact of photoelectric absorption by 
the cool unshocked winds is very small at the energies of 
the Fe xxv lines considered here. Photoelectric absorp¬ 
tion mainly affects the lines at longer wavelengths that 
are modelled by Henley et al. (2003). Its net effect is to 
skew these lines towards the blue. In our simulations, the 
only slight asymmetries between the blue and red wings 
are seen for 0 = 45 and 135°, where we note a slight de¬ 
pression in the red wing which stems from occultation of 
receding material by the star that is in front. 

Another difference with the results of Henley et al. 
(2003) concerns the fact that, in their simulations, the un¬ 
absorbed profiles of the lighter elements display a double- 
peaked morphology for viewing angles away from conjunc¬ 
tion. The reason for this behaviour is that the emissivities 
of the Lya lines of lighter elements peak at some distance 
away from the binary axis, leading to an emission region 
that has essentially the shape of a cone with its apex trun¬ 
cated (this is similar to the situation of the optical line pro¬ 
files modelled by Liihrs, 1997). Conversely, the emissivity 
of Fe xxv reaches its maximum near the apex and pro¬ 
gressively decreases along the cone. This situation leads 
to profiles that exhibit a single maximum peak near the 
line-of-sight velocity of the apex. 

We finally note that the strong line broadening for 
viewing angles away from conjunction leads to severe blend¬ 
ing of the various Fe xxv components. The resulting pro¬ 
files are thus rather complex (see Fig. 5). In practice, this 
result highlights the need to implement a least-squares de- 
convolution technique (LSD, Donati & Collier Cameron, 
1997) to recover the profiles of individual lines from ob¬ 
served data. In the remainder of this paper, we will mainly 
focus on the individual line profiles, assuming that they 
have been extracted via LSD. 

A different situation prevails for 77 = 0.3 and 77 = 0.1 
(Cases H and HI, respectively). At conjunction phases, 
we still observe a relatively narrow line, though it is now 
much wider than in Case 1. We also note that the pro¬ 
file, although narrow, is now markedly skewed towards 
the red for 0 = 0 ° and towards the blue for 0 = 180°. 
This is of course because the wind interaction zone is now 
wrapped around the star with the weaker wind (the sec- 
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Figure 4: Individual line profiles as a function of viewing angle © for the three colliding wind systems listed in Table 1. The angle 0 is related 
to the orbital inclination and the true anomaly via equation 6. 


ondary) and the post-shock material is thus moving away 
from the observer at 0 = 0° and moving towards the ob¬ 
server at 0 = 180°. Again, the maximum width (for the 
hve values of 0 considered here) is observed at quadrature 
(0 = 90°), where our calculations predict a symmetric 
profile centered on the rest wavelength of the lines. At 
intermediate phases, highly asymmetric line profiles with 
a marked peak near zero velocity are predicted, especially 
for Case III. This peak results from the fact that at these 
values of 0 and for the specific value of 0ao, one of the 
arms of the interaction region is almost perpendicular to 
the observer’s sightline. 

As a next step, we have simulated phase-resolved syn¬ 
thetic line profiles and spectra for the three cases above, 
but now assuming binary systems seen under an orbital in¬ 
clination of * = 75°. At this stage, we also account for the 
orbital motion by requesting the wind interaction zone to 


rotate with the stars on their orbit. The resulting profiles 
of individual lines at ten distinct orbital phases are shown 
in the top row of Fig. 6. The variations of the line profiles 
are clearly seen. In the following subsection, we analyse 
these profile variations by means of the so-called Doppler 
tomography technique. 

4-3. Doppler tomography 

Doppler tomography is an elegant and powerful tech¬ 
nique to study the dynamics of gas flows in spectroscopic 
binaries with circular orbits (Horne, 1991; Kaitchuck et ah, 
1994). This method is nowadays widely applied to optical 
spectra, but to the best of our knowledge, it has so-far 
never been used in the X-ray domain. 

Doppler tomography translates phase-locked emission 
line profile variations into a map of the line formation re¬ 
gion in velocity space. To do so, one adopts a reference 
frame centered on the centre of mass of the binary with 
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Figure 5: The Fe xxv line complex as a function of viewing angle 0 for the three colliding wind systems listed in Table 1. 


the a;-axis pointing from the secondary to the primary and 
the positive j/-axis pointing along the direction of the pri¬ 
mary’s orbital motion. The fundamental assumption of 
the Doppler tomography technique is that the phase de¬ 
pendence of the radial velocity u(<l>) of any cell of emitting 
gas that is stationary in the rotating frame of reference of 
the binary can be described by the simple relation: 

u(4>) = —Ux cos(2 7r$)-hUy sin(2 7r$) (7) 

Here $ stands for the orbital phase, with $ = 0 at con¬ 
junction with the primary star in front. The pair {vx,Vy) 
yields the velocity coordinates of the gas cell projected 
along the x and y axes: Vx = Vx sini and Vy = Vy sinz, 
where 14 and Vy are the actual velocities in the orbital 
plane of the binary and i is the orbital inclination. 

The Doppler map DM(vx,Vy) yields a measure of the 
flux that is carried across the line profile for each particular 
{vx, Vy) pair. If we adopt the back-projection method, the 


Doppler map can be expressed as 


DM{Vx,Vy) = 


f D(v(vx,Vy,^),^)W(^)d^ 


( 8 ) 


/ lD($)dd> 

where D(v, $) is the observed flux density at radial veloc¬ 
ity v(vx,Vy,^) (given by equation?) at phase 4. W{^) is 
the weight assigned to the observation at phase $ (see e.g. 
Horne, 1991, for a detailed discussion of the method). 

The implementation of the Doppler tomography that 
we use here is adapted from the one used by Rauw et al. 
(2005). The method is based on a Fourier filtered back- 
projection algorithm (Horne, 1991). The Point Spread 
Function (PSF) of a pure back-projection technique in the 
{vx^Vy) plane has a Gaussian core with extended wings 
having a 1 /+ Vy profile (Horne, 1991). To sharpen 
the PSF we apply a Fourier filter to the spectra prior to 
the back projection. For this purpose, we first compute 
the Fourier transform of the trailed spectrogram and we 
then multiply the result by a filter ^ exp 


to 
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Figure 6: Top row: phase-resolved line profiles for the resonance line of the Fe xxv complex for the three systems listed in Table 1 seen under 
i = 75° (from left to right: Case I, Case II and Case III). Bottom row: corresponding Doppler maps computed from 20 simulated line profiles 
for each system. The levels of the various contours correspond to 90, 80, 70, 60, 50 and 40% of the peak intensity in the map. The Roche lobes 
in velocity space are shown by the dashed contours, whilst the crosses indicate the orbital velocity amplitudes of the primary and secondary 
stars. 


suppress the 1/ y prevent the amplifica¬ 

tion of high-frequency noise in the spectra, where wjv is the 
Nyquist frequency of the spectra, ujc is set by the spec¬ 
tral resolution (for further details see e.g. Horne, 1991). 
Finally, an inverse Fourier transformation is performed to 
recover the filtered spectra that are then back-projected. 

Back-projection produces stripes across the Doppler 
map at angles corresponding to the binary phases sam¬ 
pled by the data. Therefore, the back-projection method 
can lead to the so-called ‘radial-spoke artefact’ if the data 
do not provide a uniform coverage of the orbital cycle. To 
avoid this problem, we have sampled our synthetic spectra 
at 20 equally-spaced orbital phases. 

The results of applying our method to the three sets 
of simulations are shown in Fig. 6. It becomes clear from 
these examples that the aspect of the Doppler map strongly 
depends on the shock opening angle and hence the value 
of the wind momentum ratio 77 . Doppler tomography of 
the Fe xxv line thus provides a sensitive diagnostics of the 
relative strengths of the winds in colliding wind binaries 
with circular orbits. 


Table 2: Parameters of the test models for eccentric binary systems. 


Parameter 

IV 

V 


Prim. 

Seco. 

Prim. 

Seco. 

M (10-®MQyr-i) 

12.0 

5.0 

57.0 

1.8 

Voo (kms“^) 

2400 

2500 

2860 

3200 

i?* (R©) 

16 

19 

13 

12 

M, (Mq) 

44 

50 

16 

41 

in 

62 

55 

e 

0.71 

0.896 


12 

192 

224.6 

44.6 

a (Rq) 

1727 

3300 

V 

0.43 

0.035 

X 

> 5.8 > 10.8 

> 2.4 

> 21.9 

0oo(“) 

105.1 

142.3 















































































































































































4-3. Systems with eccentric orbits 

In this section we use our code to simulate line profiles 
for a sample of eccentric colliding wind binary systems. In 
addition to the wind parameters, the orbital inclination 
and major axis, such systems are characterized by two 
new parameters: the eccentricity e and the longitude of 
periastron lu. 

Table 2 lists the parameters of the systems that we con¬ 
sider here. They are directly inspired from Cyg OB2 ^9 
(Case IV, see Naze et ah, 2012) and WR140 (Case V, see 
Pollock et ah, 2005; Fahed et ah, 2011). 

For the eccentric systems, we predict strong variations 
of the integrated line flux with orbital phase. In fact, 
since we are dealing with simulations that assume adia¬ 
batic wind interaction zones, we actually recover the 1/d 
flux variation that is expected for such systems (Stevens 
et ah, 1992). 

Figure? illustrates the results of our simulations for 
the two eccentric systems. For Case IV (e = 0.71, t] = 
0.43), the shock opening angle is relatively large, leading to 
rather broad lines at most orbital phases. The individual 
profiles have a single-peaked morphology. Their intensity, 
centroid and skewness all change with orbital phase. The 
full Fe XXV complex displays a rather complicated mor¬ 
phology indicating that LSD might again be needed to 
disentangle the variations of the individual profiles from 
apparent variations due to the blending. 

For Case V (e = 0.896, ij = 0.035), the primary wind 
overwhelms that of the secondary and, in our simple model, 
it is only the wide separation of the stars that allows for 
the shock region to remain detached from the secondary’s 
surface. The adopted pre-shock velocities are very high, 
leading to post-shock temperatures near 139 MK. This is 
above the temperature of maximum emissivity of the Fe 
XXV complex. Hence, in our simulation, the maximum of 
the emission of the line occurs at positions slightly away 
from the shock apex, thereby leading to a double-peaked 
line morphology (see Fig. 7), as discussed in Sect. 4.1. This 
renders the resulting profiles of the full Fe xxv blend even 
more complex. 

4 . 4 . Comparison with observations 

Current instruments in X-ray astrophysics lack the spec¬ 
tral resolution in the hard energy band needed to confront 
the profiles obtained in our simulations with real observa¬ 
tions. Yet, we can compare the relative variations of the 
integrated line fluxes of several eccentric systems with the 
predictions of models for adiabatic wind interaction zones. 
For this purpose, we have considered five systems which 
exhibit a clear Fe xxv feature in their broadband X-ray 
spectra. In increasing order of the orbital period, these are 
Cyg OB2 #8a (Porb = 21.9 days, e = 0.24 De Becker et ah, 
2004), WR21a (Porb = 31.7days, e = 0.69 Tramper et ah, 
2015), WR25 (Porb = 208days, e = 0.48 Gamen et ah, 
2006), CygOB2 #9 (Porb = 2.36yrs, e = 0.71 Naze et ah, 
2012), and WR140 (Porb = 7.9yrs, e = 0.89 Fahed et ah. 


2011). We have analysed XMM-Newton spectra for these 
systems. For Cyg OB2 #8a, Cyg OB2 #9 and WR 21a, the 
description of the data is given by Cazorla et al. (2014), 
Naze et al. (2012), and Cosset & Naze (2015), respec¬ 
tively. For WR140 and WR25, the data were retrieved 
from the XMM-Newton archive and processed with the 
Science Analysis Software version 14.0. The EPIC spec¬ 
tra of each object were fitted between 6 and 8keV using a 
powerlaw with zero slope for the continuum and a single 
Gaussian for the Fe K line^. The fluxes of the Fe K line are 
displayed as a function of orbital phase in Fig. 8. In this 
figure, they are compared against the l/r trend expected 
for an adiabatic wind-wind interaction zone. 

As we can see on this figure, there are clear devia¬ 
tions from the l/r trend. For the shortest period systems, 
CygOB2 #8a and WR21a, these deviations are likely due 
to the shocked gas becoming radiative around periastron 
passage and/or the shock collapsing onto the surface of the 
star with the weaker wind (Cazorla et ah, 2014; Cosset & 
Naze, 2015). The agreement is better for the longer pe¬ 
riod systems, although also here there are significant devi¬ 
ations. For instance, in the case of WR 140, the data follow 
the relation rather well up to $ = 0.986, but fall short of 
the expected emission level afterwards. WR140, CygOB2 
// 8a and CygOB2 #9 are prominent non-thermal radio 
emitters and, in the particular case of WR 140, Pittard & 
Dougherty (2006) suggested that the deviations from the 
l/r behaviour for the overall X-ray flux could reflect the 
impact of particle acceleration on the energy budget and 
the properties of the shocks (see also the case of 9 Sgr 
discussed by Rauw et ah, 2015). Whatever the exact rea¬ 
son for the behaviour seen in Fig. 8, it is clear that high- 
resolution X-ray spectra around the Fe K with next gen¬ 
eration observatories will shed new light on the properties 
of the wind interaction zone. 

5. Conclusions and future prospects 

In this paper, we have used a rather simple model, 
based on the analytical solution of Canto et al. (1996), 
to predict the morphology of the Fe xxv lines in the X- 
ray spectra of colliding wind massive binaries in the adi¬ 
abatic regime. Our results are in qualitative agreement 
with those of Henley et al. (2003). The latter authors 
computed synthetic line profiles based on the density and 
velocity fields from snapshots of 2-D axisymmetric hydro- 
dynamical simulations of adiabatic wind collisions and as¬ 
suming that the winds collide at their terminal velocity. 
Compared to Henley et al. (2003), the advantage of our 
method is its simplicity and low computational cost, al¬ 
lowing to compute large grids of models for comparison 
with actual observations. 

We have shown that the morphology of the Fe xxv 
line and its orbital changes provide direct diagnostics of 


■^The spectral resolution of the EPIC instrument is such that the 
full Fe xxv complex can be represented by a single Gaussian. 
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Figure 7: Top left: sketch of the orbit of the elliptical system of Case IV. The star with the stronger wind is depicted by the big blue dot, 
whilst the small dots illustrate the position of the secondary at six orbital phases. The shape of the contact discontinuity computed according 
to the Canto et al. (1996) model is shown at periastron and apastron. The observer is looking up from the bottom at an inclination oii = 62°. 
Top middle: synthetic individual profiles computed for Case IV and for the orbital phases shown on the left. Top right: synthetic Fe xxv 
profiles for Case IV at the same orbital phases. Bottom: same as top, but for Case V. The orbital inclination is 55° in this case. For clarity, 
all panels of the middle and right have the same vertical scale except for the periastron phase. 


the colliding wind interaction and thus of the properties 
of the stellar winds. This line offers the cleanest probe 
of the conditions near the apex of the shock region. On 
the one hand, it does not suffer from significant absorption 
by the cool unshocked winds. On the other hand, it is not 
affected by contributions from the intrinsic emission of the 
stars that make up the binary system. In this regard, there 
is a great potential for observational studies of colliding 
wind binary systems with the bolometric spectrographs 
that will fly on the coming X-ray observatories Astro-H 
and Athena. 

In the future, we will try to generalize this work also 
to systems where radiative cooling in the shocked winds 
is efficient and where the Coriolis force leads to a signif¬ 
icant aberration. An interesting case of such a system is 
V 444 Cyg (WN5 -I- 06, Porb = 4.21 days) which should be 
largely in the radiative regime and exhibits a rather promi¬ 
nent Fe xxv line in its XMM-Newton spectra (Lomax et 
ah, 2015). These properties make V444 Cyg an ideal case 
to apply the Doppler tomography technique that we have 
outlined above. 
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